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C$ Abstract. We investigate the application of nanoscale topgates on exfoliated bilayer 

*G graphene to define quantum dot devices. At temperatures below 500 mK the 

conductance underneath the grounded gates is suppressed, which we attribute to 
nearest neighbour hopping and strain. The gate-layout can thus be used to define 
resistive regions by tuning into the corresponding temperature range. We use this 
method to define a quantum dot structure in bilayer graphene showing Coulomb 
blockade oscillations consistent with the gate layout. 
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Since the first synthesis of graphene [Tj by exfoliation, the realization of quantum 
dots (QDs) in this material, e.g. as promising host for spin qubits [2j, has been an 
interesting but challenging topic [3j. One of the outstanding challenges to fabricate 
reliable graphene QDs include the difficulty in opening a sizeable and well-defined 
bandgap in graphene in order to define a QD confinement potential. It is possible 
to open a bandgap in three ways (see, e.g., [4J): (1) constrict graphene in one dimension 
to form nanoribbons [5], [6], (2) by biasing bilayer graphene [Tj, [HI E] and (3) by applying 
strain to graphene [10l EU [12] • Nanoribbon-QDs (1) produced by etching have been 
realized showing Coulomb blockade at sufficiently low temperature [13] as well as single- 
electron pumping [14J. A disadvantage of this method is that the required widths are 
lithography challenging and cause large carrier-mobility degradation. Another way to 
obtain nanoribbons utilizes unzipping carbon nanotubes which yield clean nanoribbon- 
QDs exhibiting Coulomb blockade, Kondo effect, clear excited states up to 20meV, and 
inelastic co-tunnelling [T5] , 

Following approach (2) QDs of arbitrary geometries and smooth tunable tunnel 
barriers may be achieved by local electrostatic gating of bilayers, as recently 
demonstrated by Allen et ah [16] in a dual-gate design. In their device a suspended 
bilayer membrane has been placed between a global backgate and small topgates, not 
touching the membrane. With an appropriately shaped top-gate a Coulomb-island was 
defined and corresponding conductance oscillations through this QD have been observed. 

In this work we employ on-substrate bilayer-graphene onto which we deposit a 
dielectric and metallic surface gates which are shaped to form a QD. No backgate 
has been employed. At sufficiently low temperature we observe Coulomb blockade 
oscillations, consistent with the charging energy and the location of the transport 
channel expected from the gate layout. At even lower temperature the conduction is 
exponentially suppressed over the whole gate voltage range. We attribute this behaviour 
to resistive regions forming underneath the gate metal, leading to hopping transport at 
sufficiently low temperature. As the graphene is globally covered with a dielectric and 
not in direct contact with the metal gate we propose strain as one of the possible origins 
for this behaviour. To support this hypothesis we perform simulations of the mechanical 
stress profile at the interface between top-gate and substrate. 

FigJIJa) shows a schematic of the device employed. The substrate of our device 
is undoped GaAs to increase the mobility at low temperatures due to the very low 
roughness, compared to the commonly used Si substrates with thermally oxidized 
Si02 on top [17]. In order to allow optical separation between graphene of different 
layer numbers a 200 nm thick GaAs/AlAs multilayer has been grown by molecular 
beam epitaxy onto the substrate which acts as a distributed Bragg reflector [HI [E] • 
Following the exfoliation technique [1] a bilayer of graphene has been identified by optical 
microscopy and subsequent Raman spectroscopy. For the device a flake of about 20 /im 
length and 4 /xm width has been chosen. It is contacted by pads of 10 nm thermally 
evaporated Titanium followed by 50 nm Gold. A dielectricum is fabricated [19J by first 
completely covering the whole device with an Al-layer of 2nm thickness by thermal 
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Figure 1. (a) Schematic of the QD and control device, (b) Source-drain current 
JsD versus bias voltage Usd at different temperatures, with all gates grounded. 
Characteristic of QD and control device are shown as blue and red traces, respectively. 
(c) Resistance R of the QD device versus temperature with all gates grounded. The 
solid line corresponds to a fit to R(T) = Ro exp (To/T) p with fit parameters To, Ro and 
p as described in the main text, (d) Reduced activation energy W versus temperature 
T on logarithmic scale, as explained in the main text. The solid lines indicate the slope 
corresponding to assumed p- values, as labeled. 



evaporation. After venting the evaporation chamber, the thin Al layer is fully oxidized. 
It protects the graphene during the subsequent high-energy process step in which a 
15 nm thick gate dielectric of AI2O3 is sputtered. Finally, topgates are deposited by 
thermal evaporation of 10 nm of Titanium and 40 nm of Gold. In analogy to standard 
GaAs and recent suspended bilayer graphene structures [16] the gate geometry has been 
chosen to define a QD with dimensions as shown in Fig. [jja). The continuous gates are 
labelled CI and C2, the split-gates SI and S2, respectively. For comparison, a control 
device with a 3 pm wide gate covering the bilayer has been fabricated. A voltage Usd is 
applied between the contacts labelled "source" and "drain" . 

Fig. [l](b) shows the source-drain current Tsd as function of bias voltage Usd 
at different temperatures for the QD device (blue curves). For comparison, the 
characteristic for the control device at T = 550 mK is shown by the red curve. All 
gates voltages are set to V. A non-linear characteristic with decreasing conductance is 
observed for the QD device as the temperature is lowered. The conduction is completely 
suppressed at the base temperature of T « 300 mK below a threshold voltage for onset 
of conduction at U t « ±9 mV. Note that this source-drain gap is not lifted for voltages of 
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±1 V applied simultaneously to all gates. At higher voltages gate-leakage sets in. Hence, 
even without the application of an electric field the data show an apparent transport 
gap. This gap is not observed in the wide-gate control-device [red trace in Fig. [jjb)]. 
Therefore, the conduction suppression is likely to be caused only when the surface gates 
are small and only within their proximity. Below we will discuss the assumption, that 
strain due to a mismatch in thermal expansion coefficients of gate metal and GaAs 
substrate is the origin of the conduction suppression. The local strain can be much 
stronger in smaller structures of the same thickness (see supplementary data). 

A conduction suppression only in the vicinity of the surface gates would leave a small 
island conducting in the center of the gate structure, which defines the QD. However, 
in order to allow probing the QD by transport measurement the temperature has to be 
raised to allow measurable transport through the structure. Fig.[2]shows a corresponding 
measurement at T = 550 mK. Here the conductance G is shown as function of voltage 
applied equally to gates SI and S2, Usi = Us2- A series of conduction peaks can be seen. 
We find that the application of a perpendicular magnetic field generally enhances the on- 
off ratio of these peaks (shown as blue trace), keeping the peak spacing approximately 
unchanged. For comparison, Fig. [2](b) shows the conduction-variation of the wide- 
gate control-device as function of gate voltage Uc, also measured at T = 550 mK. The 
conduction reduces smoothly towards the charge neutrality point, showing the typical 
characteristic of gated bilayer devices. No resonance or Coulomb blockade features 
have been observed. For voltages beyond 10 V gate-leakage sets in which prevented the 
observation of carrier-type change. From this plot we conclude that the transport in the 
gated graphene layer is determined by holes. 

Figs. [2](c) and (d) plot the conduction as the pairs of gate voltage (£/ci, ^02) and 
(Usi: US2) are varied, respectively. Apparently, all gates couple to the transport channel 
so that the total-conductance determining channel must be localized and close to the 
center of the gate structure. This is a strong indication that resistive regions are forming 
underneath the gate metal, confining the transport to the QD region. It can also be 
seen that while gates SI and S2 couple approximately equally to the transport channel, 
CI and C2 show a different coupling. The ratio of the corresponding gate capacitance 
is Cci/C C 2 ~ 4.2 and C S i/C S2 ~ 0.84. 

A discrete Fourier analysis of the conduction trace is shown in Fig. [21(e). It has 
been carried out for the conductance G as function of simultaneously varied Uci and 
Uq2, over a range of t/ici+C2 = — 1 . . . + 1 V. The series of power spectra was calculated 
for fixed voltages C/si+S2 simultaneously applied to gates SI and S2, over a range of 
^si+S2 = — 3 . . . + 3 V. The inset of Fig.[2](e) shows the corresponding G([/si+S2, ^01+02)- 
plot, and the main graph displays the average over all Fourier power spectra. A 
pronounced peak at a frequency of 17 V -1 can be seen, with a corresponding period 
of A[/ C i+C2 ~ 59 mV. 

Fig. [2](f ) shows the variation of the current through the device as function of Uqi 
and voltage applied between source and drain, Usb- The current variation shows a 
typical Coulomb blockade oscillation (CBO) with a charging energy extracted from 
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Figure 2. Coulomb blockade characteristic of QD device at temperature T = 550 mK. 
(a) Conductance G versus E/si+S2j for perpendicular magnetic fields of B = (red) 
and B = 3.3 T (blue), respectively. The magnetic field was set to B = 3.3 T in the 
following plots, (b) Conductance of the control- device versus gate voltage Uc- (c) 
and (d) show the conductance versus different gate voltages, as labeled. The line-cut 
at the bottom corresponds to the lowest voltage on the vertical axis, (e) displays the 
power spectrum obtained from a discrete Fourier transform of G(?7ci+C2)> averaged 
over E/si+S2 = —0.3... + 0.3V. The inset shows the conductance G as function of 
Uc\j r c2 and Z7si+S2- The upper part of (f) displays the current /sd as function of bias 
voltage Z7sd and gate voltage Uqi with all other gates grounded. The bottom part 
shows the current I$d for the cut indicated by the dashed red line in the upper panel. 



the Coulomb-diamond of Eq ~ 2meV. The resulting total capacitance Ct of the 
corresponding Coulomb island is then Ct = e 2 /Ec ~ 80 aF. From the self-capacitance 
model one would obtain a disc-shaped Coulomb island with a diameter of rfi S i an d < 
Cr/4eoe r ~ 220 nm, using a dielectric constant e r = 10.4, which is the average between 
^GaAs = 12.9 of GaAs and eAi 2 o 3 = 7.9 of AI2O3 [12]. The estimated size is consistent 
with the size of the gate defined island. Assuming that the main peak in the Fourier 
spectrum of Fig. [2](e) corresponds to Coulomb blockade oscillations of a single island the 
gate capacitance Cci+C2 = e/A[/ cl+C2 « 2.7 aF. Using the ratio Cci+C2/Csi+S2 = 2.25 
extracted from the data in the inset of Fig. [2](e) one obtaines a total gate capacitance 
of C G ^3.9aF. 

In order to investigate the nature of the interface between Coulomb island and 
electron reservoir the temperature dependence of the conduction has been analyzed 
further in Fig. [l](c) and (d). The voltage settings on the gates (all V) correspond 
to a region where conductance is not suppressed by Coulomb blockade. Therefore the 
suppression of conduction for \Usb\ < 9mV at T = 300 mK [Fig. [jjb)] originates from 
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the interface, i.e. the region underneath the surface-gates. The temperature dependence 
can be described by the exponential law R(T) = i? exp (T /T) p , where T ,R and p 
are model specific constants. In order to determine the exponent p we calculate the 
value from In W = A — plnX, where W = — <9 In i?(T)/<9 In T = p(T Q /T) p is the reduced 
activation energy and A is a constant [20]. The result can be seen in Fig.[l](d) extracting 
p — 1.1, R = 2.9xl0 6 Ohm and T = 5.3 K. An exponent p close to 1 suggests activated 
transport with an activation energy E a = 2/c B 7o = 0.9 meV, with fc B Boltzmann's 
constant. 

Considering that the conduction at T « 300 mK remains suppressed for voltages 
^ci+C2+si+S2 = — 1 V. . . + 1 V equally applied to all surface gates [see also Fig. [2](e) for 
CBO within E/ci+C2 = — IV. .. + lVatT = 550 mK] the corresponding transport gap 
must be significantly larger than the activation energy E a . Therefore we interpret the 
simple activated behavior as a signature of nearest neighbor hopping (NNH) through 
localized states within the transport gap [211 E2] . Previously, NNH has been identified in 
bilayer graphene in which a band gap has been opened by a perpendicular electric field 
in a similar temperature range [23]. For our temperature study in Fig. [l](c) and (d) no 
external perpendicular electric field has been applied that could justify the existence of a 
band gap [3 [8] . There might be a built-in electric field due to possible strong hole doping 
as indicated by the control-device in Fig. [2](b). However, this field is unlikely to lead to 
the QD definition according to the gate layout. As the graphene is globally covered with 
a dielectric and not in direct contact with the metal gate we propose strain as one of the 
possible origins for this behaviour. It has previously been established that deformation 
- described by a strain tensor s - can induce a bandgap as well [2H [25j [26]. The 
bandstructure of strained bilayer graphene has been theoretically investigated [27, 30j, 
and a band gap has been predicted for perpendicular strain exceeding s zz = 0.25. Also 
Choi et al. [31] predict a gap when the two layers of bilayer graphene are differently 
strained. Indeed, at the substrate-graphene-dielectric interface of our sample both strain 
types exist, which we conclude from a simulation of the deformation and the forces built 
up inside the material described by a stress tensor a. 

The origin of deformation and internal forces in our sample is the fact that GaAs 
substrate and the metallic gates have different thermal expansion coefficients a. The 
sample holder was kept close to room temperature during the evaporation of the surface 
gates. The temperature of the sample during the AI2O3 deposition was determined to 
be approximately 470 K. As the sample-temperature is changed the different materials 
will contract or expand at different rates and the deposited layers will excert a force 
on the substate resulting in mutual deformation. The deformation in turn generates 
an internal elastic stress a that tends to restore the material to its original undeformed 
state. 

The simulated deformation of the sample at T = 300 mK is shown in Fig. [31(a), 
with parameters provided in the caption. The simulation has been performed using 
a commercial finite element solver [28]. Plotting the displacement of the GaAs-Al 2 03 
interface (i.e. the location of the graphene flake) in Fig. |3](b) and (c) reveals that only 
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Figure 3. Simulation of sample-deformation at T = 300 mK, using thermal expansion 
coefficients a Au = 14.2 • 10~ 6 K" 1 , a Ti = 8.6 • 10~ 6 K" 1 , a A i 2 o 3 = 6.5 • 10~ 6 K" 1 , and 
c^GaAs = 5.7 • 10 _6 K _1 ; assumed substrat temperature during deposition of AI2O3 
was Tdep = 470 K; for Ti and Au Td ep = 290 K was assumed, (a) Overview of the 
deformations, scaled by a factor of 30. (b) and (c) plot the total displacement of the 
graphene layer. 



a small rectangular island bordered by CI and C2, as well as SI and S2 remains free 
of displacement and therefore stress. In contrast, the graphene-portion covered by the 
gates is strongly displaced and therefore experiences strain due to internal forces of the 
surrounding material. In addition, the bending will result in a different strain for the 
two graphene layers. 

In conclusion, we have investigated the functionality of nanoscale topgates on 
bilayer graphene exfoliated on a GaAs substrate. At low temperature the conductance 
through the gate covered graphene reduces exponentially even without the application 
of external voltages. We attributed this phenomenon to nearest neighbor hopping 
underneath the gate covered regions and propose mechanical stress imposed by the 
gates onto the graphene as a possible origin. For a certain temperature range allowing 
sufficient conduction Coulomb blockade dominated transport is observed, consistent 
with the nanoscale metallization layout. 

This work has been supported by the SPP-1459 of the German Research Foundation 
(DFG). We acknowledge the help with insulator fabrication by R. Wendisch and the 
fruitful discussions with L. Schweitzer, B. Trauzettel, B. Verberck, P. Recher, C. 
Stampfer and T. Bjorkman. 
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Appendix A: Supplementary Material 

In this supplementary material we explore the possibility of strain engineering in a 
typical top-gate structure by exploiting different thermal coefficients. The application 
of Ti-Au top-gates as used in the main text for bilayer graphene is very common in 
conventional semiconductor structures, such as for defining nanostructures in GaAs 
electron gases. The possibility of stress and induced piezoelectric fields in the latter 
type has been considered previously for instance in Ref. [29], where ideally it should be 
minimized. In the case of strain engineered graphene the opposite regime is desirable. 
In Fig. A-l[ a) we evaluate the stress tensor a in the vicinity of the gates as used in 



the main text. It shows the cross-section of a single-gate with dimensions corresponding 
to our experimental setup and the deformation being enhanced by a factor of 30. The 
corresponding stress profile of the two tensor components a zz and a xz are encoded in 
color. The location of the graphene bilayer is shown by the dashed line. For the normal 
component a zz we find tensile stress (a zz > 0, red) near the gate corner, while towards 
the center there will be compressive stress [a zz < 0, blue). In both cases strain induced 
band structure modifications and band gaps have been predicted [30j . The right part of 
Fig. A-l[a) shows the amount of shear stress a xz obtained at the position of the bilayer. 



Indeed, band gap opening by shear strain has been predicted by Choi et. al. [3]]. We 
also note that GaAs is piezoelectric, leading to electric stray fields originating from 
the strained material. Hence, the bilayer might also be exposed to an electric field 
influencing its band structure. 

Our simulation further shows that an increase in gate width W relaxes the stress, 
as illustrated in Fig. A-l[ b). It shows the stress as function of position from the gate 



center, Xcog, as its width is varied. To allow comparison, the distance was normalized 
by the corresponding gate width W. Layer thicknesses according to the measured device 
in the main text have been used in the simulations [see Fig. A-l[ a)]. It can be seen that 



at W = 350 nm the stress under the gate has relaxed to less than half of the maximum 
value at W = 130 nm of a™ ax « -80 MPa. The stress of a zz « -75 MPa in the middle of 
the gate for W = 100 nm as used in our experiment lies close to the maximum (dashed 
horizontal line). These findings are consistent with the strain-based interpretation of 
our experimental results in the main text, where the wide-gate control device remained 
conducting at low temperatures. 

The above discussion suggests that sandwiching graphene between strained 
substrate-gate structures provides a complementary tool for defining graphene- 
nanodevices with tailored transport properties. An appropriate choice of material- 
sequence, evaporation technique, thickness, geometry etc. should allow the tailoring 
of strain profile for various applications [12] . The wide opportunities of this technique 
can be derived from Fig. A-l| (d). For example, the stress can be maximized by choosing 



an optimal gold-layer thickness £a u for a given gate width of W — 100 nm, as seen 
in Fig. A-l[ c). A maximal value of a zz = —83 MPa is predicted for £a u = 27.8 nm. 



The metal thickness in our experiment was £a u = 40 nm, as shown by the horizontal 
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Figure A-l. Simulation of the deformation and stress a at T = 300 mK of a single gate 
due to different thermal expansion coefficients. Geometrical paremters are provided in 
the cross-section in (a). The deformation has been scaled up by a factor of 30. The 
colors encode the stress distributions <j zz (left part of the gate) and a xz (right part 
of the gate), as explained in the main text. The dashed orange line illustrates the 
location of graphene. (b) plots the stress felt by the graphene layer, a zz , as function 
of distance from the gate center and gate widths. The distance is normalized by the 
gate width W. (c) displays a zz as the distance from the gate center is varied on the 
horizontal axis, for fixed W = 100 nm. On the vertical axis the gold-layer thickness 
was varied. The dashed lines indicate the actual experimental conditions. In (d) <j zz 
felt by the graphene layer under the center of the gate is plotted as function of gate 
width W and gold-layer thickness. The white cross indicates the actual experimental 
condition. Between two contour lines a zz varies by 2 MPa. From the black contour 
line on the subsequent lines are blanked out because of a too close spacing. 



dashed line. The relaxation for increasing t ) however, saturates at a zz = —67 MPa for 
^Au > 80 nm. This is expected as the stress relaxes for regions in large distance from 
the substrate. 

Fig. A-l[ d) plots a zz in the center of the gate as function of £a u and W. It can 
be seen that a unique combination of W = 110 nm and <a u = 35 nm should globally 
maximize the stress. This example illustrates a possible procedure in finding optimal 
t and W parameters. The white cross in Fig. A-l[ d) indicates the stress value for the 
experimentally chosen parameters. The procedure may be extended to include arbitrary 
gate shapes, spacial thickness variations or material combinations in order to engineer 
a large variety of strain profiles. 
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